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We discuss the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive 
primordial stars. We construct a theoretical model to calculate the rate and detectability of these 
GRBs taking into account all important feedback and recent results from numerical simulations 
of pristine gas. We expect to observe a maximum of N < 0.2 GRBs per year integrated over at 
z > 6 with Swift and N < 10 GRBs per year integrated over at z > 6 with EXIST (assuming a 
sensitivity up to 10 x higher than Swift). 
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1. Introduction 



The first stars (Pop III) in the Universe are likely to have played a major role in the early 
cosmic evolution, by emitting the first light and producing the first heavy elements [fj], ^j. It is of 
great importance to understand the origin and evolution of such objects, since their detection would 
permit to probe the pristine regions of the Universe. 

Pop III stars may produce collapsar gamma-ray bursts (GRBs) whose total isotropic energy 
could be ss 2 orders of magnitude larger than average [§, @, [| . Even if a Pop III star has a supergiant 
hydrogen envelope, the GRB jet can break out of it due to the long-lasting accretion of the envelope 
itself 



2. GRB redshift distribution 

The following is a brief overview of the model we used. Readers are encouraged to check [§] 
for a full description. As long GRBs are expected to follow the death of very massive stars, their 
rate provide an useful probe for cosmic star formation history [f§ ^, 10, 11]. 



Over a particular time interval, A? t> s , in the observer rest frame, the number of observed GRBs 
originating between redshifts z and z + dz is 

, = ~j Tlbeam HGRB ^ (z) z : T~ / p{L)d\ogL, (2.1) 

dz 47T 1+z dz yiogi /im ( z ) 

where T7grb is the GRB formation efficiency, *P is the cosmic star formation rate (SFR) density, 
dV /dz is the comoving volume element per redshift unit, p(L) is the luminosity function, T[b eam is 
the beaming factor of the burst (~ 6°) and £l b s is the field of view of the experiment. The adopted 
values of Q. obs are 1.4, 2, 4, and 5 for Swift, SVOM, JANUS, and EXIST, respectively. The intrinsic 
GRB rate (the total number of GRBs per year on the sky ) is defined as = tJgrbI'(z)^^. 

2.1 Star Formation History 

In order to determine the star formation rate (SFR) at early epochs, we assume that stars are 
formed in collapsed dark matter halos. In what follows, we adopt the Sheth-Tormen mass function 
to estimate the number of dark matter halos, n ST (M,z), with mass less than M per comoving volume 
at a given redshift. The collapsed fraction of mass, F co i(z), available for Pop III star formation is 
given by 



col 



— [ dMMn sT (M,z). (2.2) 

Om JMh 2 

Where p„, is the total mass density of the background Universe. The star formation efficiency in the 
early Universe largely depends on the ability of a primordial gas to cool and condense. Hydrogen 
molecules (H2) are the primary coolant in a gas in small mass "minihalos", and can be dissociated 
by soft ultraviolet radiation. Thus, a ultraviolet background in the Lyman- Werner (LW) bands can 
easily suppress the star formation inside mini-halos. We model the dissociation effect by setting the 



lower limit of equation (2.2) as the minimum halo mass capable of cooling by molecular hydrogen 



in the presence of a Lyman-Werner (LW) background (see Hq] for more details). 
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The star formation density can be described by the following expression, 

where Q. m and are the matter and baryon cosmological density parameters respectively. The 
factors fr and Qhii represent the global filling fraction of metals via galactic winds and the volume 
filling fraction of ionized regions as shortly described below. The star formation efficiency ranges 



from /* = 0.001 as a conservative choice [ ]12[ ] to /* =0.1 [ ]13[ ] as an upper limit. Given the large 
uncertainties on the above quantities we derive upper limits on the intrinsic GRB rate using ~ 1 
year timescale radio variability surveys, which sets upper limits on the number of transients sources 
in the sky. Such limits can be used to ruled out our optimistic case. This is done by relating the 
intrinsic GRBs rate derived by our model with the number of afterglows in the sky, which is limited 
by observational limits from radio surveys (we are not really discussing how the constraints are 
derived here due to space limitation, please see [JSj] for details). Reionization: Inside growing HII 
regions, the gas is highly ionized and the temperature is ~ 10 4 K. The volume filling factor of 
ionized regions, <2hii(z), influences the typical mass of collapsed stars. Metal Enrichment : The 
metal-enrichment in the inter-galactic medium (IGM) determines when the formation of primordial 
stars is terminated (locally) and switches from the Pop III mode to a more conventional mode of 
star formation. We follow the metal-enriched wind propagation outward from a central galaxy in 
order to evaluate the average metallicity over cosmic scales as a function of redshift. The top panel 
of Fig. H shows the upper limit for Pop III SFR adopted here. 

2.2 Initial Mass Function and GRB Formation Efficiency 

The stellar initial mass function (IMF) is critically important to determine the Pop III GRB 
rate. The IMF determines the fraction of stars with minimum mass which is able to trigger GRBs. 
The GRB formation efficiency factor per stellar mass can be written as 

rM up pM up 
T7grb=/grb/ <j>(m)dm/ m<p(m)dm, (2.4) 

The /grb factor gives the fraction of stars in this range of mass that will produce GRBs. And 
(j>(m) is the stellar IMF which, we considered a Gaussian. We assume M = 4OM and dispersion 
o c = (M — Mj ow ) /3. M [ow and M up are the minimum and maximum mass adopted for Pop III stars. 
Which ranges from 1OM to 100M Q respectively. Mqrb is the minimum mass necessary to trigger 



GRBs, which we set to be 25M e [13]. 



2.3 Luminosity function 

The number of GRBs detectable by any given instrument depends on the instrument-specific 
flux sensitivity threshold and also on the intrinsic isotropic luminosity function of GRBs. We adopt 



a power-law distribution function similar to [14] 



p(L)=[£\ L<L*, (±-Y L>L*, (2.5) 
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Figure 1: Top: Optimistic model for Pop III star formation rate. The light points are independent SFR 
determinations compiled from the literature. 

Bottom: The intrinsic GRB rate (Wgrb/cIz. In other words, the number of GRBs per year on the 
sky (on-axis + off-axis). 



where is the characteristic isotropic luminosity, a = — 0.2+q 2 and j3 = —1.41^. The Pop III 
GRBs are assumed to be energetic with isotropic kinetic energy Ei so ~ \0 56 ~ 51 erg but long-lived, 
T 90 ~ 1000 s. Thus the luminosity would be moderate, L* ~ e y x 10 56 ~ 57 /1000 ~ 10 52 ~ 53 ergs/s if 
Sy ~ 0.1 is the conversion efficiency from the jet kinetic energy to gamma rays [Q]. 

Using equation ( |23| ) we can predict the observable GRB rate for the Swift, SVOM, JANUS, 
and EXIST missions. For Swift, we set a bolometric energy flux limit Fy, m = 1 .2 x 10~ 8 erg cm -2 s~ 1 . 
We adopt a similar limit for SVOM. For JANUS, Fy m ~ 10 _8 erg cm -2 s _1 . The luminosity thresh- 
old is then L y m = 47td^ Fy m . Here di is the luminosity distance for the adopted ACDM cosmology. 
Given that EXIST is expected to be ~ 7 — 10 times more sensitive than Swift, we set its sensitivity 
threshold ten times lower than Swift's as an approximate estimative. For simplicity, we assume 
that the spectral energy distribution (SED) peaks at X-to-y ray energy (detector bandwidth) as an 
optimistic case. We expect to observe a maximum of N < 10 GRBs per year integrated over at z > 



4 



Pop III GRBs forecast 



Rafael S. de Souza 



10.00 : — < — < — < — i — < — < — < — i — < — < — < — i — < — < — < — i — < — < — < — i — r 

5.00 - 




4 6 8 10 12 14 



z 

Figure 2: Predicted Pop III observed GRB rate, dA^ B /dz. Those observed by Swift, dashed red line; 
S VOM, dot-dashed black line; JANUS, dotted blue line; and EXIST, green line; for our model with Gaussian 
IMF, /* = 0.03, fcRB = 0.001. 



6 with EXIST, and N < 0.2 GRBs per year integrated over at z > 6 with Swift. 
3. Summary 

GRBs can be observed at large distances, been a powerful tool to probe the early Universe 
in situ. Together with pair instability supernovae, they are probably the most promising way to 
observe the first stars. Here we discuss the observed rate of such objects in current and future 
surveys. We take into account the recent developments concerning the origin of Pop III stars, in 



particular the last results from [15], which show Pop III stars can be less massive than previously 
thought. We expect to observe a maximum of N < 10 GRBs per year integrated over at z > 6 
with EXIST, and N < 0.2 GRBs per year integrated over at z > 6 with Swift. If we consider 
that we have no observation so far from Pop III GRBs among all GRBs observed by Swift with 
redshift measurements. We can place an upper limit over the fraction of Pop III GRBs over Pop ll/l 
ones, < l/Nswift, where Ns w ift is the number of long GRBs with redshift measurements (~ 150). 
Which gives an fraction of ~ 0.0067. Thus , considering an average of 150/6 GRBs yrs~ l (25 
GRBs yrs~ l ) with redshift measurements. It gives un upper limit of N < 0.56 Pop III GRBs yr~ l 
(0.0067 x 25 x 1/0.3), taking into account that we only measure redshifts for ~ 0.3 of all GRBs. 
Which is consistent with our estimative. 

It's important to note that even if such objects are above the flux limit threshold, most of 
the high-z GRBs might be hidden by dust and be "dark grbs". An alternative is to detect their 
afterglows. As their afterglows are expected to be long lasting ~ 1000 days, will probably be 
difficult to be detected them as transients. To overcome this problem, we should rely in surveys 
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strategies that take snapshots from the sky during long and spaced periods of time. Several tools for 
photometric supernova classification already exist and have been tested with synthetic supernovae 
light curves JT^]. We suggest that the same approach could be done in order to look for Pop III 
GRBs candidates in radio surveys. 
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